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ABSTRACT 

We study the formation of galaxies in a large volume (50ft.~^ Mpc, 2 x 288^ 
particles) cosmological simulation, evolved using the entropy and energy conserving 
smoothed particle hydrodynamics (SPH) code GADGET-2. Most of the baryonic mass 
in galaxies of all masses is originally acquired through filamentary "cold mode" ac- 
cretion of gas that was never shock heated to its halo virial temperature , confirming 
the k ey feature of our earlier results obtained with a different SPH code (jKeres et alj 
l2005f ). Atmospheres of hot, virialized gas develop in halos above 2 — 3 x IO-'^^Mq, a 
transition mass that is nearly constant from z = 3 to z = 0. Cold accretion per- 
sists in halos above the transition mass, especially at z > 2. It dominates the growth 
of galaxies in low mass halos at all times, and it is the main driver of the cosmic 
star formation history. Our results suggest that the cooling of shock-heated, virialized 
gas, which has been the focus of many analytic models of galaxy growth spanning 
more than three decades, might be a relatively minor element of galaxy formation. At 
high redshifts, satellite galaxies have gas accretion rates similar to central galaxies of 
the same baryonic mass, but at z < 1 the accretion rates of low mass satellites are 
well below those of comparable central galaxies. Relative to our earlier simulations, 
the GADGET-2simulations predict much lower rates of "hot mode" accretion from the 
virialized gas component. Hot accretion rates compete with cold accretion rates near 
the transition mass, but only at z < 1. Hot accretion is inefficient in halos above 
^ 5 X 1O^^M0, with typical rates lower than IMoyr"-'^ at z < 1, even though our 
simulation does not include AGN heating or other forms of "preventive" feedback. 
Instead, the accretion rates are low because the inner density profiles of hot gas in 
these halos are shallow, with long associated cooling times. The cooling recipes typi- 
cally used in semi-analytic models can overestimate the accretion rates in these halos 
by orders of magnitude, so these models may overemphasize the role of preventive 
feedback in producing observed galaxy masses and colors. A fraction of the massive 
halos develop cuspy profiles and significant cooling rates between z — 1 and z = 0, a 
redshift trend similar to the observed trend in the frequency of cooling flow clusters. 

Key words: cooling flows — feedback — cluster - galaxies: evolution — galaxies: 
formation — models: semi-analytic — models: numerical 



1 INTRODUCTION 

Tlie wealth of high quality observational galaxy data at dif- 
ferent epochs, from the earliest times to the present, has 
rapidly grown in the last several years. This provides a great 
opportunity to calibrate and improve our theoretical under- 
standing of galaxy formation and evolution. One of the re- 



maining challenges in theoretical extragalactic astronomy is 
the origin of the bimodality in galactic properties: massive 
galaxies that are typically red, elliptical systems with very 
little star formation, and lower mass gala xies that are typi- 
cally blue, disky, star form ing systems fe.g. lKauffmann et al.l 
I2OO3I : I Baldrv et al.ll2004l l. Usually one invokes various feed- 
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back processes to shut down star formation in massive galax- 
ies, whic h enables them to p opulate the red sequence (see 
review in lHopkins et al.ll2008l ). In addition, without feedback 
in the low mass galaxies, a large fraction of baryons would 
cool and form stars, even at very early times l|White fc ReesI 
Il978l ). contra ry to the small fraction of baryons in observed 
galaxies fe.g. iBell et al... 2003 1. 

To understand where, when, and to what degree we 
need feedback during galaxy formation and evolution, we 
first have to understand how galaxies are supplied with 
baryons, and in particular with the gas that provides the 
fuel for star formation. In our previous work l|Keres et al.l 
I2OO5I . K05 hereafter), we showed that it is the smooth accre- 
tion of intergalactic gas that dominates the global galactic 
gas supply, not accretion by merging, and that it proceeds 
via two stages. First, gas is accreted through filamentary 
streams, w here it remains re latively cold before it reaches 
the galaxy l|Katz et all 12001 K05). This process, which we 
call cold mode accretion, is very efficient, and because the 
gas does not need to cool it falls in on approximately a free- 
fall time. This means that the baryonic growth closely tracks 
the growth of the dark matter halo, albeit with a slight time 
delay. Cold mode accretion dominates the global growth of 
galaxies at high redshifts and the growth of lower mass ob- 
jects at late times. As the dark matter halo grows larger, a 
larger fraction of the infalling material shock heats to tem- 
peratures close to the virial temperature. In the denser, cen- 
tral regions a fraction of this hot gas is able to cool. This 
process, virial shock heating followed by radiative cooling, 
is the one h ighlighted in classic papers on galaxy formation 
theory (e.g.lRees fc Ostrikedl 19771 : Isii5ll977l : IWhite fc Ree3 
1 19781 : IWhite fc FrenkI Il99ll '). and to differentiate it from 
the previous mode we call it hot mode accretion. In mas- 
sive halos, this hot mode accretion not only supplies mas- 
sive galaxies with fresh g as but could also be identified as 
cluster cooling flows (e.g. lFabianlll994l ). although observa- 
tions suggest that cooling in most clusters does not actually 
reach down to galactic temperatures, perhaps because of 
heating by AGN or some other feedback mechanism (e.g. 
iMcNamara fc Nulsen|[2007l ) . The dominant accretion mode 
depends on halo mass, and the transition between these two 
regimes occurs at a mass of about Mhaio = 2 — 3 x 10^^ M©. 
At higher redshifts, cold filaments are able to survive within 
halos above this transition mass, i.e. halos dominated by hot 
halo gas. Overall, the bulk of the baryonic mass in galaxies 
is accreted through the cold accretion mode. 

Cold accretion owes its existence to the short cooling 
times present in low mass halos near the virial radius, which 
prevents the development of a stable accretion shock. The 
detailed criteria required to prevent the formation of a virial 
shock i n low mass idealized systems have be en derived by 
iBinnevI l|l977h and iBirnboim fc Dekell (120031 ). However, in 
cosmological simulations the situation is more complex than 
in these simple models. For example, the spher i cal sym me- 
try assumed in the case of iBirnboini fc Dekell (|2003D will 
not be valid owing to the presence of dense, cold filaments 
that dramatically enhance the density contrasts in the gas 
and shorten the cooling times. In halos with masses near 
the transition mass from cold to hot mode, this allows cold 
filamentary flows to still supply the central galaxy with gas 
even though some of the infalling gas shock heats to near 
the virial temperature. 



There are two main classes of galaxy formation feedback 
processes: those that prevent gas from entering a galaxy in 
the first place, "preventive" feedback, and those that ex- 
pel a fraction of the gas that does manage to enter the 
galaxy, "ejective" feedback. AGN "radio mode" heating is 
an example of pre ventive feedback (e.g. iBest et all I2OO5I : 
ICroton et al.l l200d ) and winds driven by supernovae are 
an example of e.iectiye fee dback (e.g. iDekel fc Silkl Il986l : 
ISpringel fc Hernauist|[2003l ). The effectiveness of these two 
feedback types will vary depending on the dominant accre- 
tion mode in the halo hosting the galaxy. Because cold mode 
halos have very little halo gas outside of the cold dense fil- 
aments, it is not clear if feedback processes can drastically 
affect the accretion of gas. Ejective feedback, such as galactic 
winds driven by supernovae, could still lower the masses of 
galaxies by expelling already accreted material. Such winds, 
however, might be stopped by the quasi-spherical, hot ha- 
los that surround hot mode galaxies. Conversely, preventive 
feedback, like the AGN radio mode, is likely only effective 
in hot mode halos, where it can prevent the hot, dilute gas 
that is in quasi-static equilibrium from cooling. Similar con- 
straints apply for the " quasar wind mode" of AGN feedback 
(|Di Matteo et al.llioosl 'l. which could possibly be a mixture 
of these two feedback modes. At first, a relatively small frac- 
tion of galaxy mass is ejected in the quasar wind, but later 
the energy released from the black hole accretion could in 
principle provide "preventive" feedback if the quasar was ac- 
tive within a halo that has a hot quasi-spherical atmosphere. 
In summary: ejective feedback is likely to be most effective in 
cold mode halos, and preventive feedback is likely to be effec- 
tive only in hot mode halos. An exception to this rule is "pre- 
heati ng" of inte rgalactic gas by photoioinzation ( Efstathioul 
1992i: iQuinnet a l 1996; Tliqul fc Weinberg 199^), or gravi- 
tational shocks (|Mo et al.ll2005l ). which is a form of preven- 
tive feedback that will mostly affect low mass, cold mode 
halos. 

Current cosmological simulations have insufficient res- 
olution (often by large margins) to predict these feedback 
processes from first principles, though they may incorpo- 
rate simplified recipes designed to achieve these effects. It 
is, therefore, important to gain a general understanding of 
when and where these two classes of feedback processes are 
needed to explain the observed properties of galaxies. In 
this paper, we examine the growth of galaxies by gas ac- 
cretion and the dependence of this process on galax y mass, 
halo mass, and redshift. In a companion paper (Ker es et al.l 
[2009", Paper II, hereafter), we compare the observable prop- 
erties of our simulated galaxies, such as stellar masses and 
specific star formation rates, to current data. 

In K05 we showed that the mass of the transition from 
cold to hot mode accretion does not depend much on the res- 
olution of the cosmological simulation. However, we showed 
that increasing the mass and spatial resolution substan- 
tially reduced the amount of hot mode accretion. This sug- 
gested that the simulations from K05 had too much hot 
mode accretion, perhaps the result of numerical problems 
inherent in the Smoothed Par ticle Hydrodynamics (SPH) 
formulation that we used (see ISpringel fc Hernguistl 12003 . 
SH02 hereafter, for a brief review). How much gas is able to 
cool from the hot atmosphere remained, therefore, an unan- 
swered ques tion in K05. In this paper, we use the GADGET- 
2SPH code (|Springelll20(35l ') because it does not suffer from 
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numerical overcooling in massive halos, as we demonstrate in 
Keres et al. (2009b, in preparation) and in. Keres ( 20071 ). and 
because it parallelizes more efficiently to large numbers of 
processors. The SPH algorithm implemented in GADGET- 
2, which simultaneously conserves entropy and energy, is 
a reliable method for simulations of cosmologically rele- 
vant volumes where galaxies are often only resolved with 
a small number of particles. In our previous simulations, 
this low resolution allowed a mixing of the cold and hot 
gas phases, which led to an artificial increase in the cool- 
ing rates and an overestimate of th e hot mode accretion 
IPearce et al.ll200ll : [Croft et aLlbOOll ). In our new GADGET- 
2simulations, we show that a large fraction of the massive, 
hot mode halos cannot cool their gas, owing to nearly uni- 
form density cores that develop in these halos. However, 
as discussed below ( >i3.4|l . GADGET-2simulations may still 
suffer numerical problems that distort gas accretion rates, 
albeit at lower magnitude; fully understanding these effects 
will require tests beyond those in this paper. Relative to 
K05, we also increase the simulation volume by about an 
order of magnitude, allowing much better statistics for high 
mass halos and galaxies. We use a similar but slightly better 
mass resolution. 

Here we present theoretical predictions of galaxy masses 
and accretion rates in cosmological simulations without 
strong feedback. The energy from supernovae pressurizes 
the interstellar medium, but it does not drive large scale 
outflows. We do not include any other feedback mechanisms 
that could prevent the accretion of gas or that could expel 
gas from galaxies in the mass range of interest in this pa- 
per. The only exception is that we include an extragalactic 
UV background, which affects the intergalactic medium and 
reduces gas accretion in the lowest mass halos. 

In the Paper II we confirm earlier findings that the 
amount of material that cools onto galactic component is 
too high in simulations such as ours, which causes galaxies 
at any mass, but especially at the low and high mass end 
to be more massive than their observed counterparts. There 
we argue that an efficient high redshift feedback in low mass 
galaxies (perhaps with the addition of low redshift AGN 
feedback) is the key for getting the correct galaxy masses. 

Many of the results in this paper and its companions 
(Paper II, and Keres et al. 2009b (in preparation)) first ap- 
peared as part of D. Ke res's PhD t hesis at University of 
Massachusetts, Amherst l|Keresll2007l l. 

In ij2] we describe the basic parameters and setup of 
our new simulations. In !j31we describe the smooth accretion 
of gas, both globally and for individual objects. In ^ we 
describe the properties of non-cooling and cooling halos and 
conclude in JSl 



2 SIMULATIONS 

We adopt a cold dark matter model dominated by a cosmo- 
logical constant, ACDM. In K05 we studied the properties of 
gas accretion in galaxies and their growth using an older ver- 
sion of the ACDM c osmology. Data from the first 3 years of 
the WMAP mission l|Spergel et al.ll2007l ) suggest much lower 
values of erg a nd flm (as anticipated by t he g alaxy cluster- 
ing an alyses of lvan den Bosch et al.l (|2003h and lTinker et al] 
l|2005h '). which slightly changes the time and redshift depen- 



dence of structure and galaxy growth. We therefore use the 
following cosmological parameters: flm ~ 0.26, Ha = 0.74, 



Ho/ (100 km 



Mpc 



— 0.71, and a primordial 
power spectrum index of n = 1.0. For the amplitude of 
the mass fluctuations we use as = 0.75, and for the bary- 
onic density we adopt fi^ — 0.044. All of these cosmolog- 
ical parameters are con sistent with recent measurements 
from the WMAP team (jSpergel et all 120071 ') and with var- 
ious large scale structure measurement^, except the pri- 
mordial power spectrum index which is slightly higher in 
our simulations. We model a 50.0^~^ Mpc comoving peri- 
odic cube using 288^ dark matter and 288^ gas particles, 
i.e. around 50 million particles in total. Gravitational forces 
are softened using a cubic spline kernel of comoving radius 
10h~^ kpc, approximately equivalent to a Plummer force 
softening of egrav = 7.2h~^ kpc. Using the naming scheme 
from K05, we will refer to this simulation as L50/288 later 
throughout the text. We occasionally compare this simula- 
tion to a GADGET-2simulation of a smaller volume simula- 
tion, L22/128. This simulation uses the same cosmological 
parameters and initial conditions as those that were used in 
K05. To test for resolution effects we also use Lll/64 and 
Lll/128 simulations, again with the same initial conditions 
as their analogs in K05 but simulated with GADGET-2. The 
parameters of all of these simulations are shown in Table [T] 
For the visualizations of cold accretion in i]3.3l below. we use 
a higher resolution simulation, L12. 5/288, which is described 
there. 

The initial conditions arc evolved using the SPH code 
GADGET-2 (|Springel 2005 ). The calculation of the gravita- 
tion al force is a combination of t he Particle Mesh algorithm 
fe.g. lHocknev fc Eastwoodlll98ll) for lar ge distances and the 
hiera rchical tree algorithm ( Barnes fc H ut 1986; H ernquistI 
Il987l ) for small dis tances. The smoothed particle hydrody - 
namics algorithm (|Lucvlll977l : ICingold fc Monaghan|[l977l ) 
used here is entropy a nd energy conserving, and it is based 
on the version used in lSpringel fc HernquistI (|2002| ). We use 
a modified public version of the GADGET-2. Modifications 
are made to include the cooling, the uniform UV background 
and the two-phase star formation algorithm. We briefiy de- 
scribe these below. 

We include all the rele vant cooling proce sses with pri- 
mordial abundances as in iKatz et al.r(|l996h . We do not 
include any metal enrichment or cooling processes associ- 
ated with heavy elements or molecular hydrogen. In all of 
these simulations we include a spatially uniform, extragalac- 
tic UV background that heats and ionizes the gas. The 
redshift distribution and spectrum of this background is 
slightly different than in K0 5. The backgroun d flux starts 
at z = 9 and is based on iHaardt fc Madaul (|200ll ). For 
mo re details about the ca l culati on of this UV background 
see lOppenheimer fc Davel (|2006h . We note, however, that 
smaller volume simulations of comparable resolution with 
our new UV background and with the version used in K05 
showed no noticeable differences in the evolution of the 
galaxy population above our resolution limit in the redshift 
range of interest in this paper, z = to 4. 

Once a gas particle reaches a density above the 



^ |http : //lambda.gsfc .nasa.gov/product/map/dr2/paraineters ■ cfm| 
(see the ACDM/All values) 
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Name 


L(/i"i Mpc) 


N 




'"■gas V-^^-^W / 


L50/288 


50 


2 X 288^ 





9 X 10^ 


L22/128 


22.22 


2 X 128^ 





1.1 X 10® 


Lll/64 


11.11 


2 X 64^ 


1 


1.1 X 10** 


Lll/128 


11.11 


2 X 128^ 


1 


1.3 X 10^ 



Table 1. Details of the simulations used in this paper. L is the comoving box size, TV is the total number of particles (dark+baryonic), 
2fln is the final redshift to which the simulation has been evolved and Mgas is the initial mass of gas particles in these simulations 



star forming threshold, star formation proceeds in a sub- 
resolution two-phase medium where supernova energy re- 
leased by type II SNe during star formation balances cold 
clou d formation and cloud eva poration by the hot medium 
as m iMcKee fc Ostriked (jlQTTh . Pressurization of the inter- 
stellar medium by SN feedback enables more stable gas rich 
disks, but it does not produce a large scale outflow, i.e. a 
galactic w ind. We use the same star for mation model param- 
eters as in lSpringel fc Hernguistl l|2003l ). which were selected 
to match the 2: = rel ation be tween star forma tion rate 
and gas surface density ijKennicu tt 1998; see also ISchmidtl 
1 19591 ). The code calculates the threshold density for star 
formation as the density where the mass weighted temper- 
ature of the two phase medium equals 10,000K. In practice 
this threshold density is constant in physical units during 
the simulation, and it corresponds to a hydrogen number 
density of Uh = 0.13cm~^. Each gas particle within the two- 
phase medium has an assigned star formation rate, but the 
actual conversi on from gaseous to stellar p articles proceeds 
stochastically dSpringe l fc Hcrn auistll2003h . This is similar 
to the algorithm in iKatz. (,1992i ), where each star particle 
takes half of the initial mass of a gas particle. 

To identify bound groups of cold, dense baryonic parti- 
cles and stars we use the program SKIeQ (see K05 for more 
details). Briefly, a galaxy identifled by SKID contains bound 
stars and gas with an overdensity p/pbaryon > 1000 and tem- 
perature T < 30, 000 K. Here, we slightly modify this crite- 
rion and apply a higher temperature threshold at densities 
where the two-phase medium develops. Such a modification 
is necessary to allow star forming two-phase medium par- 
ticles to be part of a SKID group, since at high densities 
the mass-weighted temperature in the two phase medium 
can be much higher than 30,000 K. To identify halos we use 
both the Friends of Friends (FOF) and Spherical Overden- 
sity (SO) algorithms with the same parameters as in K05, 
adjusted to the new cosmology. 



3 ACCRETION RATES 

The baryonic masses Afgai of SKID galaxies include 
stars, star-forming (two-phase medium) gas, and gas with 
p/Pbaryon > 1000 and tempcratur e T < 30, 000 K. Using the 
convergence tests in lMur ali ot al,J (|2002f ) as a guide, we adopt 
Afgai ~ 64mspH as our nominal resolution limit, above which 
galaxy baryonic masses are reasonably well converged. For 
the L50/288 simulation, mgpH = 9.1 x IO^Mq, making the 
resolution criterion Afgai > 5.8 x IO^Mq. To define accretion 
rates we require a galaxy to be resolved at both ends of the 

^ |http: //www-hpcc . astro . Washington. edu/tools/skld .html | 



time interval used to calculate the rate. This leaves us with 
around 10, 000 resolved galaxies at z < 2. To avoid counting 
the accretion of sub-resolution groups as a smooth accretion, 
we define smooth gas accretion as the accretion of gas parti- 
cles that were not part of any SKID identified group at the 
previous time, i.e. even bound groups smaller than 64 parti- 
cles will not be counted as smooth gas accretion. In practice 
this procedure will avoid counting sub-resolution mergers 
down to galaxies with ~ 20 — 30 particles, i.e. ~ 2 x lO^M©, 
approximately the value where the galaxy mass function in 
our L50/288 simulation drops rapidly owing to our limited 
resolution and the UV background. 

As in K05, we follow the temperature history of the 
accreted particles, noting the maximum temperature a gas 
particle reaches before it becomes part of a galaxy, Tmax- 
We ignore the temperature while a particle is in the two- 
phase medium for gas particles that are "recycled" from one 
galaxy to another. Following K05, we also split the smooth 
accretion into two components, that with Tmax < 250, OOOK, 
representing cold mode gas accretion, and that with higher 
Tmax, representing hot mode accretion. K05 found this to be 
a good empirical division between the cold and hot accretion 
modes, and they argued that a cut in physical temperature 
describes the simulations better than a cut in scaled tem- 
perature Tmax/Tvir. 

3.1 Global accretion rates 

We calculate the global smooth gas accretion rates as the 
sum of gas accreted by all resolved galaxies at a given time 
divided by the simulation volume. Following K05 we show 
results for the Tmax distribution of the global accretion in 
Figure [T] Global accretion is dominated by the cold mode 
accretion peak, i.e. by gas accreted with Tmax < 250, OOOK, 
at all redshifts including z = 0. The temperature minimum 
that separates cold from hot mode accretion washes out af- 
ter z ^ 2, caused by the very weak hot mode accretion 
in this new simulation, as we discuss below. As mentioned 
earlier, SPH codes without second-order entropy conserva- 
tion can greatly overestimate hot mode accretion (SH02). 
This problem was hinted at in the high resolution simula- 
tions evolved with P-TreeSPH in K05. A careful comparison 
between simulations that use the same cosmological param- 
eters and initial conditions but different SPH formulations 
shows that the much lower hot mode accretion rates we find 
in this new simulation are mostly the consequence of the 
different SPH formulation (Keres et al. 2009b, in prepara- 
tion). We see that the global accretion rate declines towards 
lower redshift, as we also saw in K05. The change in cos- 
mology also changes the magnitude of the global accretion 
rates, making high redshift accretion less pronounced than 
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in K05 because low as causes structure to form later. Later 
halo formation, slightly higher mass resolution, and a higher 
baryon-to-dark-matter ratio all contribute to the higher cold 
mode accretion rates seen at low redshift. 

The change in magnitude of the accretion rates with 
time is better illustrated in Figure [2j where we show the 
total smooth gas accretion rates, the cold mode accretion 
rates, and the hot mode accretion rates. The accretion rates 
peak around z = 3 followed by a drop to z = of about a 
factor of 10. The high redshift accretion rates at 2: ~ 4 — 5 
are lower by about a factor of 2-3 in L50/288 compared to 
K05. The much lower accretion rates result from the later 
formation time of halos at a given mass, which in turn results 
from the lower value of ag in the new WMAP cosmology. At 
high redshifts, this later halo formation time results in a 
smaller number of galaxies that are included in the resolved 
galaxy sample. Since high redshift accretion is completely 
dominated by efficient cold mode accretion, the magnitude 
of these rates approximately reflects the growth of the dark 
matter halos that host the resolved galaxies (minus the small 
contribution from sub-resolution mergers) , as demonstrated 
in K05. 

The fast increase in the number of galaxies that cross 
the resolution limit causes the sharp rise in global accretion 
from z~7 — 8toz~3. (Note, however, that we only count 
accretion onto pre-existing galaxies, so the first appearance 
of a galaxy above the threshold does not contribute to the 
accretion rate.) As time increases, the typical accretion rate 
at fixed galaxy mass declines, but at high redshifts this trend 
is overwhelmed by the larger number of resolved galaxies 
and by the increase in average galaxy mass (with higher 
mass galaxies having higher mean accretion rates). At lower 
redshifts, the number of resolved galaxies does not change 
significantly, and the decrease in the average accretion rate 
at fixed galaxy mass causes a drop in the global accretion 
rate. 

In Figure [5] we also compare the rate of mass growth 
through mergers to that through smooth accretion. Here, 
we include as mergers all the material that is added to re- 
solved galaxies between two simulation outputs that is not 
smoothly accreted. Almost all of the material that is counted 
as mergers joins the galaxies in dense baryonic clumps. Sim- 
ilar to the findings in K05, mergers dominate the total mass 
growth at low redshift, in this case after around z — 1. The 
gaseous part of the merger rates is now significantly higher 
than in our previous work. This is largely caused by count- 
ing the accretion of small, under-resolved, gas rich galaxies 
as mergers and in part by the later formation of halos and 
galaxies, which keeps them slightly more gas rich at late 
times. However, the gas added through mergers is always 
significantly lower than that added through smooth accre- 
tion. In a global sense, therefore, galaxies are supplied with 
gas, the fuel for star formation, mainly by smooth accretion. 

As in K05, the global star formation rate closely fol- 
lows the smooth gas accretion rate with a relatively short 
time delay. This is a consequence of the continuous gas sup- 
ply and the relatively short, observationaUy motivated star 
formation time scale we adopt here. On a galaxy by galaxy 
basis, there are systematic offsets from this tight relation be- 
tween the star formation rate and the smooth gas accretion 
rate, which we will discuss later in this section (see ij3.7p . In- 
terestingly, the z — global SFR is within a factor of two of 



the observed value (|Hopk ms fc Beacom|[200^ even though 
we do not include either supernova winds or AGN feedback. 
At higher redshifts the observational estimates have an ex- 
tremely wide range, and while our predicted star formation 
rates lie within with this range, the uncertainty of the com- 
parisons prevents strong constraints on the models. 

However, the observed star formation history is quite 
uncertain and, therefore, matching it at some given time is 
not necessarily a major success of a given model. A more 
constraining measure of previous baryonic mass accumula- 
tion is the amount of baryons that are locked in the galac- 
tic component. We plot the mass accumulation history of 
baryons in resolved galaxies as a fraction of baryonic den- 
sity of the universe in Figure [S] At z = 0, in our L50/288 
simulation around 20 percent of available baryons reside in 
galaxies and around 17 percent of all baryons are locked 
in stellar component (in and out of galaxies). For the re- 
solved galaxies only, these numbers are 19 and 16 percents, 
respectively. This is factor of 3 higher than what is observed 
(Bell et al. 2003), which suggests that, globally, too much 
gas cools onto galaxies in our simulations. This problem is 
the main topic of Paper II, where we explore this discrep- 
ancy more quantitatively and suggest what kind of feedback 
is needed to reconcile these differences. Here we only note 
that the amount of cooled baryons in the simulation ana- 
lyzed here is lower that what was suggested in earlier work 
based on SPH sim ulations without ejective feedback (e.g. 
I Benson et al.|[200ll ). A large part of this difference owes to 
the use of the entropy and energy conserving code that pre- 
vents numerical overcooling in our paper, but additional ef- 
fects likely owe to the lack of metal cooling in our simulation, 
the inclusion of a UV background which suppresses the for- 
mation of the smallest objects, and a lower as which delays 
the formation times of halos. 

3.2 Virialization of the halo gas 

Before we proceed to discuss the accretion rates of individ- 
ual galaxies in our simulation, it is interesting to check if the 
transition mass between halos dominated by cold filamen- 
tary flows and hot virialized gas occurs at a mass similar 
to that in K05, i.e. 2 - 3 x lO^M©. For this purpose. Fig- 
ure [3] plots the fraction of halo gas (excluding gas in SKID- 
identified galaxies) that is below the temperature that we 
use to divide cold and hot mode accretion, 250,000 K. We 
see that cold halo gas remains the dominant component in 
halos with masses lower than 2 — 3 x 10^^ M0, while for 
higher mass halos hot gas begins to dominate, and the most 
massive halos are almost completely filled with hot gas. The 
transition mass between cold and hot halos is almost con- 
stant with redshift even though the virial temperature of 
a halo of a given mass increases with redshift roughly as 
Tvir{z) oc (1 -|- z). One can also see that hot, virialized gas 
builds up gradually, with the transition occurring over more 
than an order of magnitude in halo mass. There is a large 
density contrast between filaments and the inter-filamentary 
region, which allows cold mode accretion to coexist with 
hot gas in the regions between the filam ents in this transi- 
ti on m ass range. The spherical model of lBirnboim fc Dekell 
(2003) predicts an approximately similar transition mass, 
smaller by a factor of 2 — 3 (though this depends on the 
precise definition of the transition in their model) . Even in 
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Figure 1. The global volume averaged smooth gas accretion rates as a function of Tmax for the L50/288 simulation. The histograms 
are calculated by adding the mass of all gas particles that were smoothly accreted in a given Tmax bin between two simulation outputs. 
Note that we express the accretion rates per decade in Tmax, unlike in K05 where they were expressed per bin. 



this idealized model, the virial shock expands outwards as 
the halo mass increases, approximately mimicking the sim- 
ulations. However, their model does not account for the fila- 
mentary structure of the infalling gas, which makes the tran- 
sition between cold and hot halos sharper than in the cos- 
mological simulations. We also compared the halo cold gas 
fractions in the simulations used in K05 with the L22/128 
and Lll/128 simulations, which used the same initial con- 
ditions but were evolved with GADGET-2, and we find very 
similar results. 

In summary, as the halo mass increases the hot viri- 
alized atmosphere slowly accumulates until it becomes the 
dominant halo gas component in massive halos. In agree- 
ment with K05, we visually confirm (see i|3.3p that cold fil- 
amentary streams survive in halos with masses much larger 
than the transition mass, especially at high redshifts, which 
makes the transition from cold to hot domination a gradual 



one and makes the col d gas fraction i n mas sive h alos slightly 
increase with redshift. lOcvirk et al.l (|2008t ) and lDekel et alj 
(2009) show similar results obtained with an adaptive re- 
finement mesh hydrodynamic code. At lower redshifts, the 
cold-to-hot transition is sharper, but it is also more sensi- 
tive to our precise definition because the virial temperature 
of the halos a t the transition mass approaches 250, 000 K. 
ICroton et all (|2006l ') suggest that the transition masses in 
K05 are biased high owing to the "geometric" SPH method 
in P-TreeSPH that was used to evolve the simulations. The 
results from the energy and entropy conserving GADGET- 
2simulation that we present here demonstrates clearly that 
this transition mass is not sensitive to the SPH method used. 
The transition mass depends, however, on the fraction of 
halo baryons available in a halo and on the gas metallicity 
as discussed by Birnboinr & Dckol (2003). In our L50/288 
simulation, in halos with IO^Mq < M;, < 5 x 10"Mq, typ- 
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Figure 2. The volume averaged smooth gas accretion, star formation and merger rates of the resolved galaxy population in the L50/288 
simulation. We show the volume averaged total smooth gas accretion rates (solid), smooth accretion in cold mode (dotted, blue) and 
smooth accretion in hot mode (dashed, red). We also plot the star formation rate densities of the resolved galaxy population (green, 
long-dashed) , the volume averaged rates of mass increase through mergers (dot-dashed) and the gaseous part of the merger rates (tripple- 
dot-dashed). The upper panel plots the rates versus cosmic time and the lower panel plots the rates versus redshift. In both panels we 
also overplot the fraction of baryons locked in resolved galaxies (triangles) and the fraction of baryons locked in stars in resolved galaxies 
(stars), with the corresponding values on the right hand side y-axis. 



ically around 40 percent of the halo baryons are residing out 
of the galactic component. Our preliminary results from sim- 
ulations that include the metal line cooling and supernova 
driven galactic winds, which increase the amount of gas in 
halos, indicate that the transition between halos dominated 
by cold and hot gas shifts upward by a factor ~ 3 in halo 
mass. 

3.3 Illustrations of Cold Mode Accretion 

How do filaments that are driving the cold mode accretion 
connect to galaxies and halos? How do these connections de- 
pend on halo mass and redshift? The answers to these ques- 
tions depend somewhat on numerical resolution. In higher 
resolution simulations, the central cores of filaments are able 
to reach higher densities and therefore cool faster and be- 
come thinner. Furthermore, the low mass galaxies and sub- 
halos that can be tracked at higher resolution can drag or 
disturb cold gas filaments in their parent halos. To investi- 



gate the qualitative behavior of filamentary cold accretion, 
we therefore use the highest resolution simulation currently 
available to us, evolved to z — 1 with 288'^ particles in 
a 12.5/i"^ Mpc comoving box (L12. 5/288). This simulation 
has no ejective feedback, but it was evolved with metal line 
cooling and stellar mass loss that gets recycled into the sur- 
rounding gas. The detailed properties of the galaxies that 
form in this simulation will be presented in future work, but 
it is worth mentioning that this extra physics does not make 
a significant difi^erence to the properties of the intergalactic 
halo gas, since without winds the metals and recycled ma- 
terial stay within or close to galaxies. We explicitly checked 
whether the buildup of the hot gas in halos changes with 
this additional physics and found no difference, i.e. the tran- 
sition mass between cold gas and hot gas dominated halos is 
the same as in all the other simulations used in this paper. 
The gravitational softening in this simulation is 220/i~^ pc 
(Plummer equivalent) physical at z = 3 and it increases to 
440/i~^ pc physical at z = 1. Such high resolution is nec- 
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Figure 3. The fraction of cold halo gas in a halo of a given mass at four different redshifts. Lines show the running median of the cold 
gas fraction (blue, solid) and the hot gas fraction (red, dashed). (We exclude any galactic gas.) 



essary to see the variety of structures in the cold accretion 
flows and to properly resolve the regions very close to the 
galactic disks. 

Figure 13] shows the characteristic structure of cold mode 
gas in halos below the transition mass: a halo at z = 3 
with Mh = 1.1 X IO^'^Mq (upper panels) and a halo of the 
same mass at 2: = 1 (lower panels). The particles are color 
coded by temperature, and their relative density on the plot 
indicates the projected column density (because the masses 
of the gas particles out of galaxies are usually identical). To 
illustrate the filamentary geometry, we show two different 
projections of each halo. 

As emphasized in K05, owing to the presence of the 
dense filamentary streams, a typical high redshift halo has 
regions of great density contrast. High densities and short 
cooling times in these fila ments make conditions unfavor- 
able for shock propagation (|Binnevll 19771 : iBirnboim fc Dekell 
I2OO3I ). In the inter-filamentary regions, where the gas den- 
sity is often one or two orders of magnitude lower, some 
gas shock heats even when the halo mass is much lower 
than the transition mass. This is in contrast with simple ID 
models, which necessarily average away density variation at 
a given radius, and which typically show a quick transition 
from cold mode to hot mode. In the upper panels of Fig- 
ure m dense filamentary fiows connect to the outskirts of 
the central galaxy and to an infalling satellite galaxy (see 



§ 13. 5|) . Tidal interaction has created a bridge between these 
two galaxies. Near i?vir the filaments are thicker than the 
galaxies themselves, but they are then compressed by the 
surrounding shock heated (virial temperature) gas, which 
extends to ~ 0.5i?vir. Near the galaxies, the filaments are 
narrower than the galactic disks. 

The situation changes at low redshift, where lower mass 
galaxies typically reside in larger filamentary structures 
comparable in cross section to the size of the halo. However, 
the densest parts of the filaments can still survive within 
the virial radius. Deeper in the halo, some of this gas gets 
shock heated to the virial temperature and some stays cold, 
resulting in a mixture of cold and hot accretion that sup- 
plies the gaseous galactic disk visible in the lower panels of 
Figure 2] 

Figure [S] shows more massive halos, with mass ~ 
IO^'^Mq, a factor ^ 3 above the transition mass seen in Fig- 
ure[3l The limited volume of our L12. 5/288 simulation does 
not contain any halo of this mass at 2 = 3 so we show halos 
at 2 = 2 and 2 = 1. Here the situation changes dramatically. 
At 2 = 2, a typical ~ W^^Mq halo evolves at the intersec- 
tion of several filamentary structures. The shock heated gas 
in this halo fills up the cavities between the filaments, and 
in some directions it extends to several virial radii. One can 
clearly see how gas that falls onto a halo from the region 
between the filaments gets heated at a shock front, typically 
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Figure 4. Temperature of the gas in low mass halos at z = 3 and z = I. Panels show the gas in a region of 250/i~^ kpc (comoving) 
on a side and 250h~^ kpc (comoving) in projected depth. The virial radius is shown as the circles. Upper panels show the two different 
projections of a halo at 2 = 3 dominated by a cold halo gas with = 1.1 X IO^^Mq, while lower panels show a halo at z = 1 with 
Mfi = 1.12 X IO^^Mq. Vectors attached to particles show their projected velocities. Note that warmer gas temperature within the densest 
clumps is caused by sub-resolution implementation of the two-phase star forming medium. 



slightly outside i?vir. The middle and the lower panels show 
only the gas with T < lO^K, to reveal the penetration and 
filamentary structure of the cold accretion. These filaments 
are enhanced by infalling substructures, which drags parts 
of the filaments into the more massive halos and often splits 
larger filamentary structure into multiple fiows separated by 
hot gas. Furthermore, halos of this mass are dominated by 
hot gas at T ~ Tvir, which provides a high pressure envi- 
ronment that compresses the filaments. This results in the 
thin filamentary streams, most of which can penetrate to the 
halo's central regions. The zoom-in panel of the z = 2 halo 
shows that these cold filaments survive to supply both the 
central and the satellite galaxies (in the upper right corner 
of this panel) with cold gas. However, a large fraction of the 
infalling gas is now shock heated to the virial temperature 
and fills the halo, so the accretion rate of the central galaxy 
in this halo is a factor of several below the halo accretion 



rate. The filamentary streams in this panel show clumps 
and disturbances from substructure. The filaments dragged 
by the infalling substructure typically infall with a larger 
impact parameter bringing in large amounts of high angular 
momentum gas. 

Results for the z — 1 halo (right hand panel) are similar, 
but the incoming filaments are somewhat thicker compared 
to J?vir because a 10^^ M0 halo is no longer an unusually high 
density peak. At fixed overdensity, the physical density and 
cooling rates are lower, so a smaller fraction of cold filament 
gas survives its journey to the halo center, and more is shock 
heated to join the hot atmosphere. 

In the most massive halo in our L12. 5/288 simulation, 
with a mass of Mh « 9 x IO^^Mq, the situation is different 
yet again as we show in Figure (6] In this high sigma peak, 
there are again multiple filaments entering the virial radius, 
and the inter-filamentary space is filled with shock heated 
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Figure 5. Temperature of the gas in Milky Way mass halos at 2 = 2 and z = 1. Panels show the gas in a region of Ih ^ Mpc (comoving) 
on a side and lh~^ Mpc (comoving) in projected depth. The virial radius is shown as the circles. Left panels show an Mh = l.lx lO^^Mg 
at 2 = 2, while the right panels show an = 1.2 X IO^'^Mq halo at z = 1. Upper panels show all gas particles. The middle panels 
show only the gas with T < K but with the same color scale (indicated on the right). The lower panels also show only the low 
temperature gas, but are zoomed-in to show a region 375/i~^ kpc (comoving) on a side (approximately within 0.8-Rvir) a-^d 12bh~^ kpc 
(comoving) deep. Vectors attached to particles show their projected velocities. Note that warmer gas temperature within the densest 
clumps in middle and lower panels is caused by sub-resolution implementation of the two-phase star forming medium. 
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Figure 6. Temperature of the gas in a group size halo at z = 1, with = 9 X IO^^Mq. Panels show the gas in a region 2/1-1 
(comoving) on a side and 2h~^ Mpc (comoving) in projected depth. The virial radius is shown as the circles. Upper panels show all gas 
particles. The middle panels show only the gas with T < 10^ K but with the same color scale (indicated on the right). The lower panels 
also show only the low temperature gas, but they are zoomed-in to show a region 7b(ih~^ kpc (comoving) on a side (approximately 
within 0.8-Rvir) and 2b0h~^ kpc (comoving) deep. Vectors attached to particles show their projected velocities. 
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gas that extends well beyond -Rvir- Compressed by this hot 
medium, the filaments are even thinner compared to the 
virial radius. However, the middle and lower panels show 
that these cold flows do not penetrate deeper than about 
0.5i?vir- Some of the filaments are still connected to satellite 
galaxies, but the inner satellites and the central galaxy are 
separated from this gas supply. The bottom, zoom-in panel 
shows that the remnants of the filamentary flows clump into 
numerous dense cold clouds. Some of these dense concentra- 
tions of gas are actual galaxies, but most are a consequence 
of fragmentation of the filamentary flow induced by shocks 
and high pressure, cooling instabilities of the warm/hot gas, 
or clouds that form out of stripped galactic material. Accre- 
tion of such clumps can provide fresh gas for central galaxies 
in massive halos, at least in simulations without strong pre- 
ventive feedback. 



3.4 Accretion rates of individual galaxies 

The accretion rates of individual galaxies as a function of 
their mass and of the central galaxies as a function of their 
parent halo meiss are shown in Figure [7] We also plot the 
median accretion rates of the central galaxies as a solid line 
and the median accretion rates of satellite galaxies at a given 
mass as a dashed line (left panels). The larger simulation 
volume, compared to K05, allows us to highlight some details 
that were previously not very clear, especially for the high 
mass end at low redshifts, which we will discuss shortly. 

At 2 = 4 all the galaxies are completely cold mode dom- 
inated, and the accretion rates have a steep dependence on 
galaxy and halo mass. A typical IO^^Mq galaxy and a typical 
central galaxy in a lO^^M© halo accrete at 10 — 15MQ/yr, 
while galaxies in halos of around 5 x lO^'^M© accrete at 
30 - lOOMo/yr.H At z; = 2 galaxies at all masses are still 
cold mode dominated, but the dependence of the accretion 
rate on galaxy and halo mass becomes very shallow. For 
a change in halo mass of two orders of magnitude, from 
IO^^Mq to IO^^Mq, the smooth gas accretion rate of cen- 
tral galaxies increases by only a factor of ~ 2. The massive 
halos at 2 = 2 are already full of hot virialized gas, and only 
a tiny fraction of this halo gas cools, resulting in the very 
low hot mode accretion rates. Between z = 2 and z = 1, the 
median accretion rates at the low mass end drop by factors 
of ~ 2 — 3. Not only are the trends with mass now shallow, 
but above ~ 2 x 10^'^ Mq the accretion rates of the central 
galaxies actually start to decrease with mass, though they 
turn up again above IO^^Mq. 

In the mass range where the accretion rates are de- 
creasing with mass, the halos are full of hot virialized gas 
that cannot cool. At the highest mass end, where the gas 
accretion rates once again increase with mass, cold mode 
accretion dominates. In massive halos only the remnants of 
the cold filaments survive at low redshift, in the form of cold 
clouds. Some of these clouds are bound and identified with 

^ The typic al accretion rates in Figure 7 are lower than those 
reported by iDekel et al. I f2009') by a factor of several. However, 
they report total gas fluxes through to the inner 15 kpc of their 
simulated halos rather than the smooth gas accretion rates of their 
central galaxies; the difference between these quantities could be 
significant. 



SKID, so we do not count them as smooth accretion even 
though they are not galaxies in the classical sense, i.e. they 
are not embedded in dark matter halos. Cold clouds not 
identified by SKID as galaxies also contribute to this cold 
mode accretion in massive halos. In addition, some cold, 
galactic gas gets decoupled from the stellar component ow- 
ing to strong ram pressure and tidal forces. 

The accretion of cold clouds with these various origins 
causes the cold mode accretion in very massive halos at late 
times. While part of this accretion could be realistic, much 
of it is probably a consequence of numerical artifacts, for 
the reasons discussed in i]4.3l below. Since most of this ac- 
cretion comes from clouds with masses slightly below the 
resolution limit (see i]3.7|l . we refer to it as "cold drizzle". 
Whether or not this cold drizzle is physical or numerical, 
the fact remains that group and cluster mass halos, which 
have masses at least an order of magnitude above the tran- 
sition mass between cold and hot mode accretion dominated 
galaxies, effectively stop accreting gas from their hot virial- 
ized halos, i.e. the hot gas no longer cools. 

The results a.t z = Q reveal some interesting trends. The 
dependence of the accretion rates on galaxy and halo mass 
are similar to z; = 1, but the rates are lower overall. At the 
low mass end the rates increase from around 0.2 — O.3M0/yr 
in halos with masses of ~ IO^^Mq to 1 — 1.5M0/yr in 
~ IO^^Mq halos. Once again, the rates drop for masses 
above lO^'^M©. At halo masses larger than several times 
10^^ M©, the halos are completely dominated by hot mode 
accretion, except for the most massive halos where poten- 
tially spurious cold drizzle provides several M0/yr of ac- 
creted gas and where the hot mode rates are even lower. 
Most halos have lower accretion rates at z = than at 
z = 1, and smooth accretion is almost completely shut off in 
halos with masses larger than several times IO^^Mq. How- 
ever, there is a sizeable population of halos with masses of 
~ 10^^ — 5 X 10^^ M0 whose accretion rates have increased 
since z = 1, and the mass dependence of the upper envelope 
of accretion rates is quite steep. These halos h ave develope d 
what could be considered classic cooling flows l|Fabianlll994l ). 
In the same mass range, the majority of the halos have ac- 
cretion rates, especially in the hot mode, that are very small. 
We will discuss the detailed properties and statistics of cool- 
ing and non-cooling halos in SJ?) 

The left panels of Figure [7] also list the fraction of galax- 
ies with extremely low accretion rates, less than 15% of the 
median accretion rate of galaxies at 2 x lO^^M©. Specifically, 
we adopt thresholds of 0.1, 0.3, 1 and 3M0/yr at z = 0, 1, 
2, and 4, respectively. These galaxies typically accrete only 
one or zero gas particles between the two simulation out- 
puts, and most such objects are satellite galaxies in larger 
parent halos. 

3.5 Gas accretion in satellite galaxies 

Figure[7]includes a comparison of central and satellite accre- 
tion rates in our GADGET-2, L50/288 simulation. At z = 4, 
the satellites and central galaxies have very similar accre- 
tion rates at a given galaxy mass. At z = 2, the satellites 
already accrete less than the central galaxies, but only a 
small percentage of galaxies have their accretion rates dras- 
tically reduced. The relative difference between central and 
satellite galaxies increases with decreasing redshift and is 
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Figure 7. Smooth gas accretion rates of individual galaxies in the L50/288 simulation. In the left panels we show the accretion rates of 
all resolved galaxies as a function of their mass and in the right panels we show the accretion rates of the central (i.e. the most massive) 
galajdes in a given halo as a function of their parent halo mass. The satellite galaxies are marked with crosses, while the central galaxies 
are plotted as points in the left panels and diamonds in the right panels. In both columns, black symbols represent total accretion rates. 
In the right panels we also show cold mode accretion in blue and hot mode accretion in red. The solid line shows the median accretion 
rate of central galaxies at a given mass while the dashed line is for satellites. In the upper right corner of the left panels we also indicate 
the fraction of resolved galaxies with extremely low accretion rates (see the text for details) .Note the different vertical scale for the z = 
panels. 
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quite large for lower mass galaxies at 2; = 1. At 2; = 0, the 
diflterences are substantial. The median accretion rate of the 
lowest mass satellite galaxies, those around the resolution 
limit, is zero. Note that owing to finite particle mass and the 
time interval between the two simulated outputs, we cannot 
measure accretion rates below 0.07Mq /yr. However, even at 
2 = around 45 percent of satellite galaxies are still able 
to accrete at rates above O.lM0/yr. The differences in the 
accretion rates of higher mass satellites and central galaxies 
of equal mass decrease with galaxy mass, but they are still 
noticeable even in satellites with masses of ~ 2-3x10^^Mq. 

At any redshift, roughly 80% of the galaxies with ex- 
tremely low accretion rates are satellites, except at z = 4 
where there is a roughly equal mix of satellite and central 
galaxies (but at this redshift the low accretion rate galaxies 
are only 1% of the total). The fraction of resolved galaxies 
with low accretion rates grows from 0.06 at 2 = 2 to 0.12 at 
2 = 1 to 0.24 at 2 = 0. Although the median accretion rate 
of all galaxies is dropping over this interval, the increase in 
extremely low accretion rate galaxies is driven entirely by 
the growing gap between the central and satellite popula- 
tions. This growing gap is also evident in the median rates 
at A/„,i 2 X lO^^MtT), the ratio of median accretion 



■ e.g., 



2 at 



1 and a factor of ~ 10 at 



rates is a factor of 

2 = a 

ISimha et al.l (|2008l ) investigate the relative accretion 
and merger rates of satellite and central galaxies in the 
L22/128, P-TreeSPH simulation that is the main simula- 
tion analyzed by K05, and they discuss implications for the 
properties of central and satellite populations as a function 
of halo mass. Although the form of our analysis is differ- 
ent, the general trends here are consist ent with thos e foun d 
in the more detailed investigation by ISimha et al.l l|2008h . 
In particular, despite the difference in simulation code, we 
find that a significant fraction of satellite galaxies have 
continuing gas accretion, and (in an investigation not il- 
lustrated here) that the satellites whose accretion is most 
strongly suppressed ten d to be low mass systems in high 
mass halos. In addition, iKeresI (|2007| ) analyzed satellite ac- 
cretion in the cosmologic al simulations using P-TreeSPH, 
GADGET-2and Gasoline l|Wadslev et all |2004 ) SPH codes 
and also found that the fraction of accreting satellites in- 
creases with increasing redshift and that the fraction of 
non-accreting satellites increases with increasing halo mass. 
These results are consistent with the view that satellites 
within the virial boundary of a halo may continue to grow 
by accreting gas from massive substructures in which they 
reside. These predictions (and the implied continuation of 
star formation in s atellite galaxies) are su pported by the 
empirical study of of lWeinmann et all (|2006l ). who find that 
satellite galaxies in the SDSS are less red than in semi- 
analytic models (SAMs) of galaxy formation. Especially at 
high redshifts, our simulation results indicate that the as- 
sumption commonly used in SAMs, that satellite galaxies 



(adopted by, e.g., .Kauffmann et al.lll993 
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I2OO3I. and in sub- 



sequent models from these groups), is not valid. At low red- 
shifts this approximation might be correct for the bulk of 
galaxies, but there seems to be a mass dependence on the 
effect of satellite strangulation, and a substantial fraction of 
satellites experience continuing accretion. 
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Figure 8. The fraction of total smooth gas accretion in cold mode 
(blue points) as a function of galaxy (left panels) and halo mass 
(right panels) at z = A, 2, 1 and 0. The lines represent the median 
cold mode (blue, solid) and hot mode (red, dashed) fractions in 
a bin 0.2 dex in mass. 



3.6 The fraction of accretion in cold and hot mode 

The relative importance of hot and cold mode accretion is 
shown in Figure |8] where we plot the fraction of the to- 
tal smooth gas accretion in these two modes as a func- 
tion of galaxy and halo mass. At 2 > 2 the accretion is 
dominated by the cold accretion mode at all masses. At 
2 = 1, cold mode dominates at the low mass end, i.e in 

-,11 
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galaxies with baryonic masses lower than 2 — 3 x 10^ 
or halo masses lower than 2 — 3 x IO^^Mq. From Afhaio ~ 
2 X IO^^Mq — IO^^Mq, cold and hot mode accretion are com- 
parable on average. Above Mhaio ~ 2 x IO'^^Mq cold mode 
once again dominates, but the spread in the cold accretion 
fraction of individual objects is very large. At 2 = 0, cold 
mode dominates at low masses, the hot mode begins to dom- 
inate at halo masses > 2 — 3 x 10^^ Mq, and it continues to 
dominate until masses of 10^^ Mq. Above this mass, the me- 
dian cold and hot mode fractions are comparable, but the 
spread from halo to halo is very large. The significant cold 
mode fraction at high halo masses arises from a combination 
of very low hot mode accretion rates and "cold drizzle" (see 
i]3.4|l . which may be mostly a numerical artifact. The abso- 
lute rates of cold accretion in these halos are still fairly small. 
We see that when a transition between cold and hot accre- 
tion occurs, the transition masses are Afgai ~ 2 — 3 x IO^^Mq 
and Afhaio ~ 2 — 3 x IO^'^Mq, very similar to the values in 
K05. However, the much lower hot accretion rates in the 
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GADGET-2simulation changes the appearance of this plot 
relative to K05's figures 5 and 6: cold mode dominates at all 
galaxy and halo masses at high redshift, and at low redshift 
the hot mode domination at intermediate masses gives way 
to (possibly spurious) cold mode domination at the highest 
masses. 

3.7 Star formation rates 

In Figure |9] we show the star formation rates of the individ- 
ual galaxies (on the left) and of the central galaxies in halos 
(on the right) as a function of the galaxy and halo mass. 
We plot the median star formation rate as a solid line and 
the median smooth accretion rate from Figure [7] as a dashed 
line. We also show the smooth gas accretion rates when we 
allow sub-resolution clumps to be counted as accretion (dot- 
ted) and the total gas accretion rates, which includes the gas 
supplied by mergers (dot-dashed). In addition we also show 
median of the total baryonic mass supply (i.e. including both 
gas and stars) . In all four cases we exclude galaxies with zero 
star formation rate. The overall star formation rate follows 
the smooth gas accretion rate, but there are systematic de- 
viations from this trend. 

In low mass galaxies at low redshift, the typical gas ac- 
cretion rates are higher than the star formation rates. The 
accreted gas in low mass galaxies needs to build up a suf- 
ficient reservoir before a significant fraction of the gas can 
cross the star forming threshold density. This physical effect 
is modulated by a numerical artifact: the gas densities in the 
central parts of these marginally resolved objects are under- 
estimated, which delays the onset of efficient star formation 
until even higher masses. Because of this gas accumulation 
at the low masses, lower mass galaxies are always much more 
gas rich than their higher mass counterparts. As many low 
mass gal axies are indeed obse rved to be more g as rich at 
both low (|BeU fc de JondbOOd ') and high redshift (|Erb et all 
I2OO6I ). part of this effect could indeed be physical, but the 
exact delay in star formation also depends on the resolution. 
The simulations with higher mass resolution show that this 
behavior shifts to lower masses as higher densities can be 
resolved, allowing the star formation rate to follow the gas 
accretion rate more closely. 

Conversely, the star formation rates of high mass galax- 
ies are typically higher than their accretion rates because of 
the larger contribution of minor and major mergers, which 
provide galaxies with gas as well as stars. In addition, the 
same objects usually have slightly higher gas accretion rates 
at earlier times, so the delay between accretion and star for- 
mation boosts the latter relative to the former. Comparison 
of the dashed and dotted lines also indicates that a large 
amount of the gas accretion in massive galaxies and halos 
comes from sub-resolution clumps (cold drizzle), which (as 
discussed previously) may be a numerical artifact. Below 
Mgai ~ 1O"M0 or Mhaio ~ lO^^Mo, adding sub-resolution 
mergers to smooth accretion makes a negligible difference. 
The gas supply from resolved mergers appears to be only 
mildly important and only at the high mass end. Mergers 
do dominate the total mass buildup of massive galaxies (e.g. 
iMaller et al.| [2006') at late times but most of the mass they 
supply is in form of stars. 

At all times galaxies form a well defined star forming 
sequence with mass, with a slope similar to the observed 



slope (|Daddi et al.ll2007l : iNoeske et all 120071 '). However, the 
normalizatio n of the sim ulated relation is lower than that 
observed fsee lDavj|2008l ). Similar to the gas accretion rates, 
the typical star formation rate for central galaxies in Milky 
Way size halos (~ lO^^M©) is ~ 2M0/yr at z = and 
grows to ~ 4OM0 /yr at z — 4, albeit with a sizeable scatter. 
The star forming sequence has less scatter as a function of 
galaxy mass than as a function of halo mass. This indicates 
that the galactic star formation law also plays a role in the 
regulation of this sequence, together with the main driver 
of gas supply, cold mode accretion. As we discuss in Paper 
II, reproducing the observed distribution of star formation 
rates in high mass galaxies requires that the gas accretion for 
most of these systems be sharply reduced but that a fraction 
of them retain the accretion and star formation predicted by 
the simulation. 

3.8 Galaxy buildup 

An intermediate or high mass galaxy may be dominated 
by hot accretion at z — 0, but it spent earlier phases of 
its growth below the cold-to-hot transition mass, and at 
high redshift where all galaxies are dominated by cold mode 
accretion. In addition, a massive galaxy may have grown 
partly by mergers with systems that were themselves built 
by cold mode accretion. Figure [TU] shows the contribution of 
gas initially accreted through cold mode to the final masses 
of galaxies at z = 0. Here, we select all the particles that are 
contained in resolved galaxies at z = and trace their tem- 
perature history back through time. For each particle, we 
find the maximum temperature that it has ever reached be- 
fore that redshift. We exclude any times when the particle's 
density was higher than the two-phase medium threshold 
in GADGET-2, to avoid counting temperatures that were 
increased through the sub-resolution star formation algo- 
rithm. If the temperature of the particle has never exceeded 
250, 000 K, we count it as being accreted through cold mode. 

As expected, a typical low mass galaxy forms com- 
pletely through cold accretion. The contribution of hot mode 
accretion increases until Mgai ~ 4 — 5 x lO^^M©, where the 
cold mode contribution is ~ 70 — 75%, but then the trend 
reverses and the hot mode contribution starts decreasing in 
more massive galaxies. This reversal is caused by the in- 
creased contribution of low mass mergers and mergers at 
early times, which both add mostly cold mode accreted ma- 
terial to the massive galaxies. The near absence of hot mode 
accretion in massive halos reinforces this trend. The median 
cold mode contribution at masses Mgai ~ 5 x lO^'^M© is 
higher than 90%, and it varies between 75% and 95% at 
even higher masses. The dispersion in the cold mode frac- 
tion is large at all masses. Removing "cold drizzle" does not 
alter our conclusion that even the most massive galaxies are 
built from baryons that were originally acquired via cold 
mode accretion. 

The increase of the cold mode contribution at the high- 
est masses is a consequence of "downsizing" in galaxy forma- 
tion, whereby progenitors of the most massive galaxies today 
formed most of their stars and acquired most of their mass at 
earlier times than less massive galaxies. This trend is estab- 
lished by numerous observations, and while it is sometimes 
described as "anti-hierarchical," it is predicted both by semi- 
analytic models and by numerical simulations based on the 
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Figure 9. Star formation rates of individual galaxies in the L50/288 simulation. In the left panels we show the star formation rates of 
all resolved galaxies as a function of their mass and in the right panels we show the star formation rates of the central (i.e. the most 
massive) galaxies in a given halo as a function of the parent halo mass. The solid line shows the median star formation rates and the 
dashed line over plots the median smooth gas accretion rates similar to Figure [7] The dotted line shows the smooth gas accretion rate if 
we include sub-resolution clumps (see text for details). The dot-dashed line includes all channels of gas supply, smooth accretion, clumpy 
accretion and gas delivered through resolved mergers, while the long-dashed line also includes the supply of stellar component (mostly 
through mergers). The galaxies with zero star formation rates are excluded. 
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Figure 10. The fraction of galactic baryonic mass initially accreted through cold mode plotted as a function of galaxy mass at 2: = 0. 
The solid line plots the median fraction acquired through cold mode in bins 0.15 dex in mass. 



ACDM cosmology. We will discuss galaxy buildup in more 
detail in Paper II. However, it is clear from our results here 
that that owing to the high efficiency of cold mode accre- 
tion in lower mass objects, the low efficiency of gas cooling in 
massive halos, and the hierarchical buildup of the high mass 
objects, the typical galaxy at any mass assembles primarily 
through cold mode accretion. While important in some in- 
dividual cases and in a limited mass range, cooling from the 
hot viriali zed atmosphere, as emphasized in the sta ndard 
paradigm l|Rees fc Qstriketll 19771 : 1 White fc ReesHlOTSl ). rep- 
resents only a secondary channel of galaxy growth. While 
feedback and metal cooling could increase cooling from the 
hot atmosphere, as discussed in the next section, these ef- 
fects would also increase the mass range over which cold 
mode accretion dominates, as discussed in H3.2I We there- 
fore expect this conclusion to hold even in simulations that 
include these processes. 



3.9 Resolution effects 

We briefly discuss the effects of resolution on our findings. 
For this purpose we use the Lll/64 and Lll/128 simula- 
tions (see Table [T|, which have the same initial conditions 



as the corresponding simulations in K05 but were evolved 
using the GADGET-2code. The mass resolution of Lll/64 is 
comparable to our L50/288 simulation, while Lll/128 has a 
factor of eight higher mass resolution (a factor of two in spa- 
tial resolution) , including additional small scale power in the 
initial conditions. We compare these simulations only down 
to 2: = 1, since structure on scales larger than the simulation 
box becomes nonlinear at lower redshifts. 

Comparing the galaxy mass function at masses that 
should be resolved by both simulations, at the low mass end 
the Lll/128 simulation has a higher mass function while at 
the high mass end they both agree quite well. The smooth 
gas accretion rates in Lll/128 relative to L64 are around 
30 percent lower for most halo masses, over the whole red- 
shift range. The same applies to the star formation rate. 
These lower rates owe to the higher resolution, which en- 
ables lower mass objects to form at earlier times, convert- 
ing part of their gas supply into stars and hence lowering 
their gas accretion rates at a given redshift. With higher 
resolution the star formation in any galaxy also begins ear- 
lier because higher gas densities can now be reached, and it 
therefore locks more baryons into stars at a given time, again 
decreasing the amount of gas available for star formation at 
later times. 
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Owing to the small volume, there are only two halos 
more massive than 3 x lO^^M© at 2; = 1 in the Lll sim- 
ulations. In the case of central galaxies in these two halos, 
differences in the star formation rate are quite small at z > 2 
but increase to 30-40 percent by z = 1. The cooling of hot 
gas in massive halos is very inefficient in both simulations, 
but it is even less efficient in LII/I28, which is not surprising 
as a slightly larger fraction of baryons are now locked up in 
galaxies. The effect of "cold drizzle" from the sub-resolution 
objects, which dominates the supply of gas to massive galax- 
ies at late times, is similar at high resolution if we count 
only accretion below the 64 x mspH,Lii/i2s threshold, but 
it is significantly less pronounced in the higher resolution 
simulation if we count the gas accretion from groups below 
64 X mspH^LU/fn- This result remains the same even if we 
account for the gas contribution from resolved mergers (as 
some of the cold clouds are now above the resolution thresh- 
old) , implying that there is a significant decrease in the total 
accretion of cold gaseous clouds in the higher resolution sim- 
ulations. At least part of this effect must owe to the fact that 
galax;ies of a given mass are actually more gas poor in the 
higher resolution simulations, and there is therefore less gas 
available for stripping once they fall into a massive halo, but 
it could also suggest that the evolution of the cold clumps is 
greatly affected by our limited resolution. We also find that 
the typical mass of the cold gas clumps in the halo are lower 
in the higher resolution simulation. 

The Lll/64 vs. Lll/128 comparison suggests that our 
conclusions about the declining importance of hot mode 
accretion in massive halos would if anything be stronger 
at higher resolution. They also indicate the robustness of 
our qualitative conclusions about the feedback mechanisms 
needed to reproduce observed galaxy masses and specific 
star formation rates, which is the subject of Paper II. The 
significant resolution dependence of the accretion and star 
formation rates for the most massive galaxies implies that 
the simulation results for these galaxies should be inter- 
preted with caution. In particular, these tests provide some 
evidence that continuing cold accretion in the most massive 
galaxies at late times, is partially a numerical artifact, a 
point we return to in i]4. 31 below. 



4 COOLING AND NON-COOLING HALOS 

In this paper, as in K05, we find that galaxy growth in ha- 
los below Mhaio ~ 3 X 10^^ Mq is dominated by cold accre- 
tion, and that dominant hot gas components develop only 
above this mass. The important difference relative to K05 is 
that the hot mode accretion rates in these massive halos are 
much lower. The sign of this difference is as expected given 
the difference in SPH algorithms used, with the entropy- 
conserving algorithm providing a better separation between 
the cold and hot gas phases. However, as we show below, 
the cooling rates in most of our simulated massive halos are 
much lower than predicted by typical semi-analytic calcula- 
tions in the absence of feedback. Furthermore, the range of 
cooling rates is large at fixed halo mass. In this section we 
discuss the origin of this large dispersion of cooling rates, 
the nature of the continuing cold accretion in high mass ha- 
los, and the possible connections to the observed properties 
of X-ray groups and massive galaxies. 



4.1 Gas profiles: cores and cusps 

The right panels of Figure [7] show the accretion rates of cen- 
tral galaxies, and here we focus on the behavior at high halo 
masses. Above 10^^ M©, there are 63 halos at z; = and 28 
halos at z = 1; at 2; = 2 we lower the threshold slightly 
to 8 X 10^^ Mq to improve statistics, yielding 3 halos. At 
2 > 1 there is no clear demarcation between cooling and 
non-cooling halos, but if we take an arbitrary division at 
3.5M0 yr~^, then 100, 65, and 13 percent of the massive ha- 
los are above this threshold at z = 2, 1, and 0, respectively. 
If we consider only hot accretion, these percentages change 
to 33, 0, and 3. 

Figure [TT] compares the density and temperature pro- 
files of two massive (Mhaio ~ 4 x lO^^M©) halos at z = 0, 
one in which the central galaxy has a high hot mode accre- 
tion rate of 3OM0 yr~^ and one with a (more typical) hot 
mode accretion rate of O.3M0 yr~^. We only show particles 
that are part of a halo but are not part of the individual 
galaxies within a halo. The density profiles (shown by plot- 
ting the densities of individual SPH particles vs. radius) are 
strikingly different: the cooling halo has a roughly power- 
law profile in to 0.01-Rvir, while the non-cooling halo has a 
core (more accurately, a sharp break to a shallower power 
law) inside ~ 0.15i?vir. The temperature of the hot gas is 
roughly constant in to 0.01i?vir in the cooling halo, but in 
the non-cooling halo it rises by about a factor of two from 
-Rvir to 0.02i?vir, iuside of which essentially all gas particles 
are cold and belong to the central galaxy. These density and 
temperature profiles are typical of other high mass halos in 
these two regimes. 

The right panels of Figure [TT] show the instantaneous 
cooling times of gas particles in the two halos. In the cool- 
ing halo, particles out to 0.1i?vir have typical cooling times 
below 10^" years. In the non-cooling halo, only a small frac- 
tion of particles have cooling times below 10^" years. We 
can predict hot mode accretion rates from these plots by 
taking the amount of gas that can cool during the typical 
time between simulation outputs (^ 1.3 Gyr at 2; ~ 0) and 
dividing by the time interval. We obtain 36M0 yr~^ and 
O.I7M0 yr~^, in comparison to the measured hot accretion 
rates of 3OM0yr~^ and O.35M0 yr~^. Thus, the difference 
in hot mode accretion rates can be understood as a direct 
consequence of the difference in density and temperature 
profiles. While dynamical processes like shocks, infall, dy- 
namical friction and turbulent motions likely contributed 
to the formation of cored profiles (e.g. Conrov fc Ostriker| 
l2008l : iDekel fc BirnboimI I2OO8I : iKhochfar fc Ostriker , 200 
and they may have injected heat that offset cooling in the 
earlier stages of halo formation, at late times the simple 
"instantaneous" estimate of cooling can explain the very 
low accretion rates found in 33] without the need for any 
additional heating processes. We briefly checked the higher 
redshift results and found that similar cores start to form 
in the most massive halos in our simulation even at 2 = 2, 
albeit with slightly higher central densities. We will explore 
the profiles of the higher redshift sample in future work. 

The density and temperature profiles of these two halos 
are reminiscent of the gas profiles in observed cooling flow 
and non-cooling flow cluster s (e.g. iDe Grandi fc Molendil 
l2002l : iMorandi fc Ettoril I2OO7I ') . This leaves the question of 
what creates the range of profiles in the first place. In the 
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Figure 11. The gas properties at 2 = of two group size halos with masses of 4 X lO^^M©, one that is not cooling as is typical (upper 
panels) and a rarer one that can cool (lower panels). The left panels plot the densities of gas in the halo (points). Model profiles are 
shown for an isothermal gas density (solid line) and the density calculated assuming an isothermal gas in hydrostatic equilibrium within 
an NFW potential for a IO^'^Mq halo, where dotted line shows the initial density profile and dashed line shows the gas density profile 
after it cools for 2 Gyr. The middle panels plot the gas temperatures and the right panels plot the gas cooling times. The right panels 
also indicate hot mode accretion rates from the simulation, an accretion rate from a simple cooling calculation, and the cooling rates 
from the SAM used in Croton et al. (2006) (see text for details). 



cluster context, a common hypothesis is that mergers dis- 
rupt the quiescent gas evolution that would otherwise pro- 
duce cooling flows. Large scale shocks during major merg- 
ers can raise the central entropy of the gas, thereby low- 
ering the central density and increasing cooling times. As 
a rough test of this hypothesis, we compared the mean as- 
sembly redshifts (defined as the redshift at which the stellar 
mass of the largest progenitor exceeded half the mass of the 
z = galaxy) of the central galaxies of massive z = halos 
with hot mode accretion rates above and below 3.5M0 yr~^. 
These assembly redshifts indicate the time at which the cen- 
tral parts of the halo were largely in place. While the scat- 
ter is large and the number of systems rather small, we find 
that for Mhaio > 8 x W^^Mq the mean assembly redshift 
of cooling halos is z ~ 1.7, while for non-cooling halos it 
is slightly below z = 1. Thus, it appears that halos with 



a long quiescent phase after the assembly of their central 
regions are most able to build cu spy density proflle s that 
allow substantial cooling (see also iBurns et al.ll2008l ). This 
can be understood from the buildup of the dark matter ha- 
los. The formation of dark matter halos proceeds in two 
stages, a rapid growth by major mergers at early times an d 
slow growth by minor mergers at late times (|Li et al.ll2007l ). 
The time when the rapid growth stage finishes increases 
with increasing halo mass. While a large fraction of halos in 
the range of 10^^ — lO^^M© typically still experience major 
mergers around z = 1, most such halos have finished their 
rapid growth hy z — 0. For more massive halos, however, a 
signiflcant fraction are still experiencing major mergers at 
< z < 1, though others are already in a stage of more 
quiescent growth. 

In this scenario, cooling flows set in when the growth of 
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halos slows. Observationally, this view is supported by the 
study of Vikhhnin et aLi 12003) > who find that none of the 
~ 20 clusters they observe in X-rays at z > 0.5 show any sig- 
nature of strong cooling flow proflles. At low redshift, how- 
ever, roughly half of the clusters in their sample have cooling 
flow proflles. In the bulk of clusters in their higher redshift 
sample, they also find some signatures of recent merger ac- 
tivity, and they conclude that the higher frequency of major 
mergers at higher redshifts is responsible for the difference 
in cooling flow fractions. In our simulation, we find no high 
mass halos with hot accretion rates above 3.5Mq yr~^ at 
z = 1. Cold accretion would not produce X-ray emission, so 
it is not relevant to this comparison. At z = 0, 3 percent 
of the M > 10^'^ Mq halos have hot accretion rates above 
3.5MQyr~^, and 17 percent above IMQyr"^. The fraction 
of cooling flows is thus low compared to the observations, 
but the redshift trend is correct. The predictions are likely to 
depend on halo mass as well as redshift, so quantitative com- 
parison to data will require larger volume simulations with 
more halos in the mass range (~ 10" - 10^^ M0) probed by 
X-ray cluster observations. 

Our preliminary tests of simulat ions with stellar wind 
feedba ck and metal cooling from lOppenheimer fc Pavel 
(2006) show that even with the larger fraction of gas in these 
halos, gas density cores are still common in a mass range sim- 
ilar to that in our simulation. The central parts of these cores 
are able to cool at slightly higher rates than in our no-wind, 
no metal simulations. Despite the low cooling rates of our 
present simulations, additional feedback in group mass halos 
might be needed to keep their massive central galaxies red 
enough to match the observations. However, the formation 
of a gas density core that prevents or at least signiflcantly 
slows down cooling from the hot atmosphere appears to be a 
generic result of hierarchical structure formation, even with- 
out feedback from AGN or other astrophysical sources. The 
exact fraction of halos affected, with or without additional 
feedback mechanisms, remains an open question. 



4.2 Implications for Semi- Analytic Modeling 

Returning to Figure 1111 the solid blue lines show singular 
isothermal sphere (SIS) gas proflles, p oc r~^, scaled to the 
virial mass and gas fraction of these halos. This density pro- 
flle describes the cooling halo fairly well, but it drastically 
overestimates the density inside Q.li?vir for the non-cooling 
halo. The solid green curve shows the density proflle for 
isothermal gas in hydrostatic equi librium in the proflle of 
an NFW halo (H avarro et al.lll997l) with t he concentration 
expected for IO^Mm (|Bullock et aLlboOll ). again scaled to 
the virial mass and gas fraction of our simulated halos. We 
modify the NFW potential to include the gravitational force 
softening at small radii, and this modification produces the 
flattening of the density proflle inside ~ 0.05-Rvir. In Keres et 
al. (2009b, in preparation), we evolve spherically symmetric 
gas distributions with these initial profiles using GADGET- 
2. The dashed green curve shows the density and tempera- 
ture profile after the isothermal NFW distribution has been 
evolved (with cooling) for 2 Gyr. The evolved density proflle 
approximately follows the SIS curve, while the temperature 

proflle drops inside ~ 0.03i?vir. 

The piecewise cooling solution of lBertschinge3 (Il989l ) 
provides accurate estimates of the gas accretion rates at halo 



center for proflles such as these l|Keresll2007f ). Semi-analytic 
models frequently use an approximation to this solution, as- 
suming either SIS profiles or isothermal gas in equilibrium in 
an NFW halo, to calculate accretion rates in massive halos. 
Comparing the cooling time profiles in Figure [TTl it is clear 
that this approach will be reasonably accurate for the cool- 
ing halo but will overestimate the accretion rate in the non- 
cooling halo drastically. As a sp ecific example , we a pplied 
the cooling algorithm used in the lCroton et al] (|2006l ) semi- 
analytic model: assuming an isothermal gas distribution and 
the gas fraction of the simulated halos, we calculate the ac- 
cretion rate as drricooi/dt = O.Smgasrcooi Vc/-Rvir) where rcooi 
is defined as tcooi(7'cooi) = -Rvir/Vc. (We ignore the small dif- 
ference in virial radius deflnition.) The predicted accretion 
rates are 45M0 yr~^ for the cooling halo, which is accurate 
at the 50% level, and 38M0 yr"'^ for the non-cooling halo, 
which is more than two orders of magnitude above the ac- 
tual accretion rate in the simulation. More detailed compar- 
ison of the simulation results to a full range of semi-analytic 
models will be presented elsewhere (Lu et al. 2008). 

It is possible that the cores found here in most high mass 
halos are enhanced by some numerical artifact caused by the 
effective "phase separation" between hot and cold gas that 
arises in the entropy-conserving formulation of SPH. How- 
ever, our tests on idealized cooling flow models (^Kere§ 2Q0j) 
suggest that GADGET-2gives reasonably accurate results for 
halo masses and particle numbers like those considered here 
(more accurate than the codes without the entropy conserv- 
ing SPH formulation). Taking the numerical results at face 
value, there are significant implications for semi-analytic 
models, since with standard assumptions they would greatly 
overestimate the cooling rates for the majority of high mass 
halos in this simulation. Some models reset their cooling ra - 
dius after a major merger (e.g. lSomerville fc Primaci3ll999l ). 
which goes in the right direction, but an isothermal post- 
merger gas proflle still leads to an overestimate of the cooling 
rates, by an amount that depends on the time interval and 
the algorithm used for the cooling calculation. Models that 
assume that th ere is a core in the gas density profile of mas- 
sive halos (e.g. iBower et al.ll2006D could in principle provide 
a better match to the central parts of our simulated halos, 
but the central temperatures in such models should also be 
increased to match the simulated gas properties. Moreover, 
models that assume a density core in all halos fail to rep- 
resent the cooling flow systems that develop at late times. 
If semi-analytic models are consistently overestimating gas 
cooling rates in massive halos, then they are likely requir- 
ing too much feedback in these halos to compensate. More 
accurate calibration of the cooling recipes in semi-analytic 
models can be developed by comparing the evolution of gas 
proflles in hydrodynamic simulations and the semi-analytic 
calculations; the short analysis here provides a flrst step in 
this direction. 



4.3 Cold Drizzle 

Even at high halo masses and low redshifts, accretion onto 
central galaxies is usually dominated by cold gas in sub- 
resolution clumps, the form of accretion that we have de- 
scribed as "cold drizzle." Thermal in stability can form cold 
clumps in situ from hot hal o gas (|Mo fc Miralda-Escudd 
1 19961 : iMaller fc BuUockl |2004) . but the cold drizzle gas in 
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our simulations has never been hot, so thermal instability 
clouds are not its source. (Our simulation probably lacks 
the resolution needed to model thermal instability in galaxy 
mass halos in any case.) Rather, the cold drizzle appears to 
be a mix of gas that has been stripped from larger galaxies 
and sub-resolution clumps that form within cold/ warm fila- 
ments that are disrupted at the outskirts of high mass halos, 
where they may be compressed by the surrounding hot gas. 

Both of these sources should produce cold gas that can 
accrete onto galaxies in massive halos. However, there are 
several numerical effects that could artificially amplify this 
form of accretion in our simulation, perhaps by large fac- 
tors. Firs t, ram pressure strip ping of cold gas from infalling 
galaxies (jGunn fc GottI Il972l ) is almost certainly overesti- 
mated at the resolution of L50/288 because the typical 
smoothing lengths of the hot particles greatly exceed the 
size of the galaxie s themselves (see tests and discussion by 
iTittlev et al.llioOll '). Second, a similar effect could cause cold 
gas clumps (regardless of origin) to sink to the center of 
the halo more quickly than they should, as artificial ram 
pressure slows their tangential velocities. Finally, SPH codes 
have difficulty following surface instabilities that should de- 
velop wh en cold dense clum ps move through a hot dilute 
medium l|Agertz et al.ll200'if ). so clouds that should be de- 
stroyed by instabilities or by conduction (which is not in- 
cluded in our simulations) may artificially survive their jour- 
neys to the halo center. As discussed in i l3.9l we find that the 
amount of cold drizzle in massive halos at z = 1 is reduced 
when the mass resolution is increased by a factor of eight rel- 
ative to that of L50/288, though the limited redshift range 
and box size of our resolution tests does not allow us to draw 
clear conclusions about the influence of resolution on cold 
accretion rates in massive halos. In a systematic comparison 
of simulatio ns run with different codes from the same initial 
conditions l|Keresll2007l ). we find more cold drizzle in mas- 
sive halos in simulations run with GADGET-2compared to 
those run with PTreeSPH and Gasoline. This difference sug- 
gests that the effects mentioned above may be more severe 
for the GADGET-2SPH formulation because of the different 
calculation of the hydrodynamic forces and the different av- 
eraging of the properties of gas on the interface between the 
dense cold phase and the hot dilute medium, which could 
result in different rates of cloud survival and formation and 
different ram pressure forces at moderate resolution. 

Our findings on low accreti on rates in massiv e halos are 
qualitatively similar to those of lNaab et al.l (|2007l ), based on 
much higher resolution GADGET-2simulations of the forma- 
tion of individual galaxies in a cosmologica l conte xt. Like the 
simulation analyzed here, the lNaab et al.l |2003) simulation 
incorporates a simple star formation algorithm and has no 
strong feedback in the form of supernova -driven wind s or 
AGN heating. At their highest resolution, iNaab et al.l find 
that the central galaxy in a halo of several times W^^Mq 
stops forming stars and is not supplied with new gas after a 
redshift z ~ 1; thus they, too, find that accretion from the 
hot gas halo shuts off without the assistance of feedback. In 
their lower resolution simulations, however, accretion con- 
tinues at lower redshift, perhaps in the form of cold drizzle 
like that in our simulations. Increased resolution could sup- 
press this cold accretion both by more accurately treating 
ram pressure and instabilities and by depleting the gas sup- 



ply in low mass clumps (which get converted into stars when 
the resolution is high enough). 

We remain uncertain whether the cold accretion rates 
in massive halos in our simulation are approximately correct 
or substantially overestimated because of finite resolution 
and numerical effects in the GADGET-2SPH formulation. 
We show in Paper II that cold drizzle makes only a small 
contribution to the final masses of galaxies in these high 
mass halos, but it can significantly affect their star formation 
rates and hence their colors. 

4.4 A Natural End to Smooth Accretion? 



Following the re sults of iKatz et al.l (|2003l ) and 
iBirnb oim fc Dekell (|2003l ) on cold and hot accretion, 
several authors noted that simulations or semi-analytic 
models of galaxy formation would better reproduce obser- 
vations of the galaxy luminosity function and the low star 
formation r ates of massiv e galax ies if hot accretion were 
suppressed (iBinnevI l200i : K05; iDekel fc BirnboimI l2006l : 
ICattaneo et al.l 2003). Motivated in part by observations 
of galaxy clusters, these authors and others suggested 
that preventive feedback from AGN might preferentially 
suppress hot accretion, since relative to cold accretion the 
hot gas has lower density, higher temperature, and larger 
geometrical cr oss-section. In o ne con crete implementation 
of this idea, ICattaneo et al.l l|2006l ) and ICattaneo et al.l 
(2008) show that a semi-analytic model that assumes the 
shut-down of accretion and star formation in halos above 
a critical mass Mhaio ~ IO'^^Mq, together with traditional 
semi-analytic model assumptions about stellar feedback, 
reproduces many of the observed features of the galaxy 
red sequence and bimodality of the galaxy popul ation. 
The semi-analytic models of ICroton et al.l l|2006l ) and 
iBower et al.l l|2006l ) allow AGN feedback to suppress gas 
accretion in halos that, according to the cooling radius 
criteri a in these models, hav e developed quasi-static hot gas 
halos (|White fc Frenklll99ll ). and these models also achieve 
better agreement with observations than earlier generations 
of semi-analytic models. 

In this paper we find that the hot accretion rates in 
simulated massive halos are low even without any strong 
feedback mechanisms, principally because of the formation 
of cored gas density profiles that have long central cooling 
times. If we assume that cold drizzle in massive halos is 
mostly a numerical artifact that would disappear at high 
resolution ( i]4.3|l . we conclude that the physics of hierarchi- 
cal halo assembly and galaxy formation on its o wn leads to 
a shu tdown of gas accretion in massive halos. iNaab et al.l 
(2007*) emphasize that the simulated galaxies that form in 
their highest resolution simulations have many of the prop- 
erties of observed massive ellipticals — in particular, no 
ongoing gas accretion and star formation. Maybe smooth 
accretion comes to a natural end, with no need for AGN 
feedback or other suppression mechanisms. 

However, the story does not appear to be quite so sim- 
ple. In our simulations, the fraction of cooling halos increases 
towards lower masses, and while we find some density cores 
even m halos of ~ 10^^ Mq, core formation is not enough to 
shut down gas cooling in most halos in the 10^^ Mq — 10^'^ M© 
range. If we remove hot mode accretion and cold drizzle in 
massive halos "by hand," then we find that the most mas- 
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sive galaxies are predominantly "red a nd dead," in agree- 
ment with the high resolution results of iNaab et alj (|200'if ). 
However, in Paper II we show that the simulated massive 
galaxies are too massive by a factor of 2 — 5 compared to 
observed galaxies of the same space density, so in the absence 
of feed back the processes simulated here and bv lNaab et al.l 
l|2007fl do not appear able to reproduce the galaxy baryonic 
mass function; at a given halo mass, galaxies in the real uni- 
verse are less massive. On the other hand, we also show in 
Paper II that reproducing observed specific star formation 
rates requires that central galaxies in a sizeable fraction of 
10^^ - 10^^ M0 halos remain on the star forming sequence, 
so a complete shutdown of accretion above IO^^Mq is not 
what is desired. Finally, we note again that ejective, galac- 
tic wind feedback in lower mass galaxies will ultimately have 
significant impact on the most massive galaxies by reducing 
the stellar and gas masses of galaxies that merge with them, 
by increasing the overall IGM gas supply, and by enriching 
the IGM and therefore enhancing cooling rates. Our findings 
here are important steps towards understanding the physics 
of galaxy growth, and they differ in some significant ways 
from conventional wisdom, but they do not provide a full 
solution to the puzzles of observed galaxy evolution. 



5 CONCLUSIONS 

Our principal results are based on a cosmological SPH sim- 
ulation of a 50^~^ Mpc (comoving) box in a ACDM uni- 
verse, evolved with 288^ dark matter and 288^ gas parti- 
cles using the code GADGET-2. The combination of volume 
and particle number allows us to resolve (with > 64 par- 
ticles) galaxies of baryonic mass Mgai > 6 X IO^Mq while 
following the formation of several cluster-size halos with 
A^haio ^ 1O^'*M0. Relative to the main simulation analyzed 
by K05, our present simulation has an eleven times larger 
volume , and it uses the entropy-co nserving SPH formula- 
tion of ISpringel fc Hernguistl (|2002| ). which should provide 
a more accurate treatment of multi-phase gas than the 
PTreeSPH code employed by K05. We incorporate photoion- 
ization by the UV background and thermal feedback on a 
2-phase galaxy interstellar medium, but we do not explicitly 
add galactic winds or track metal enrichment, so this is a 
"minimal feedback" simulation. 

Confirming the key result of K05, we find that galaxies 
with baryonic masses below 2 — 3 x 10^" or halo masses 
below 2-3 X 10" Mq gain most of their mass through "cold" 
accretion of gas that h as never been close to the halo virial 
temperature (se e also iKav et ah I I2OOOI : iFardal et al.l I2001I : 
iKatz et all 120031 '). In halos below ~ 2.5 x IO^Mq, the ma- 
jority of intergalactic gas is colder than 250,000 K, while 
in higher mass halos the majority of gas is close to the 
virial temperature; the transition mass stays nearly con- 
stant from 2: = 3 to z = 0. Simi lar behavior is found in 
spherically sym metric calculations l|Birnboim fc Dekeill2003l : 
iDekel fc Birnboim 2006 ). and the likely cause of the transi- 
tion is that the short coohng times in low mass halos prevent 
the f ormation of a stable vi r ial shock ([B irnboim fc Dekell 
l2003l : see also iBinnevI [l977l : IWhite fc F^ cnk 1991). How- 
ever, as shown by K05 and illustrated more comprehensively 
here, halos near and above the transition mass have hot 
gas halos penetrated by cold gas filaments, which may con- 



tinue to feed the central galaxy. Similar results are found in 
adaptive mesh re finement simulations (lOcvirk et ahl I2OO8I : 
iDekel et al.ll2009l : see also iDekel fc BirnboimI (|2006l ) figure 
6). At high redshift these cold filaments penetrate close to 
the halo center, while in high mass halos at 2; < 1 they 
are usually truncated or disrupted before reaching the cen- 
ter. The difference likely reflects the higher physical densi- 
ties and shorter associated cooling times of the high red- 
shift fllaments that prevent the expansion of a stable shock 
(Birnboim fc Dckcl 2003) in these high density regions. 

The key difference between our present results and those 
of K05 is that hot accretion rates are much lower, a conse- 
quence of changing to the more accurate, entropy-conserving 
SPH formulation. Because of the low hot accretion rates, 
the growth of high redshift galaxies [z > 2) is dominated by 
cold accretion at all galaxy and halo masses. At lower red- 
shifts, massive halos experience a combination of hot accre- 
tion and "cold drizzle," accreting lumps of cold gas below our 
64-particle galaxy resolution thresho ld. Our resolution tests 
here and those of iNaab et al.l (|2007l ) suggest that some or 
even most of this "cold drizzle" may be a numerical artifact. 
Eliminating it would significantly affect the predicted star 
formation rates and colors of massive galaxies, but it would 
not have much impact on their predicted stellar masses (Pa- 
per II). At all galaxy masses, the majority of galaxy baryonic 
mass in our simulation was originally acquired by cold ac- 
cretion, either directly onto the main progenitor or onto a 
satellite that merged with it. 

Galaxy accretion rates decline systematically from z = 
4 to z = 2 to z = 1 to 2 = 0. For 10^°Mq < Mg^i < lO^M©, 
there is a trend of increasing accretion rate with increasing 
mass at all redshifts, though it becomes flatter at low z. At 
higher masses, accretion rates have large scatter and are sen- 
sitive to cold drizzle. With our adopted star formation pre- 
scrip tion, empirically motivated by the Kennicutt-Schmidt 
law (iKennicuttll 19981 : ISchmidtll 19591 ). star formation closely 
tracks gas accretion with a short delay; K05 showed this on 
a global basis, and here we show it on a galaxy-by-galaxy 
basis. At each redshift, more massive galaxies tend to have 
higher star formation rates but (moderately) lower specific 
star formation rates. Satellite galaxies orbiting within the 
virial radius of larger halos experience continuing accretion 
and star formation. The differences in the accretion rates 
of central and satellite galaxies of the same mass are small 
at high redshift, while at low redshift a significant fraction 
of satellites have much lower accretion rates. More detailed 
analysis of the satellite accretion and merger rates ( based 
on the K05 simulation) appears in ISimha et al.1 (|2008i ). 

The low hot accretion rates in our simulated massive 
halos arise because these halos typically have cores in their 
hot gas density profiles (or, more precisely, because they 
break from an approximately profile to a much shal- 
lower profile in their inner regions). By z = 0, a modest 
fraction of the M > IO'^^Mq halos develop cuspy density 
profiles with cooling cores, evolution that is reminiscent of 
Vikhlinin et al.'s (2007) finding that cooling flow clusters 
are fairly common at z = but much rarer at 2 > 0.5. 
We speculate that the cores are produced by shocks dur- 
ing the chaotic, rapid assembly phase of halos, and that the 
more quiescent evolution during the slow accretion phase 
allows cooling to establish a cusped profile. The accretion 
rates of central galaxies in our massive halos are predicted 
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with reasonable accuracy by instantaneous cooling calcu- 
lations that use the actual simulated gas density profiles, 
but cooling calculations that assume isothermal gas in equi- 
librium in an NFW halo or singular isothermal halo may 
overestimate the cooling rates by orders of magnitude. Our 
simulation predicts low hot accretion rates even though it 
does not include AGN heating or other forms of preventive 
feedback. While such feedback may be needed to explain the 
absence of cold gas in coo ling flow clusters (i Kaastra et al.l 
I2OOII : [Peterson et al.ll2003 ). semi-analytic models may over- 
estimate its global importance by assuming incorrect halo 
gas profiles (and thus excessive cooling). 

The comparison to observations in Paper II shows that 
this simulation, despite its low hot accretion rates, over- 
predicts the galaxy baryonic mass function. The amount of 
baryons that is locked in galaxies at 2: = and the amount 
of baryons locked in stars at 2: = is about a factor of 2-3 
higher than observed. This is caused by simulated galax- 
ies being too massive at all masses, but the exact differ- 
ence is a function of galaxy mass with the lowest and the 
highest mass galaxies being the most problematic. Match- 
ing observations therefore requires either the suppression of 
cold accretion (which seems unlikely on physical grounds) or 
substantial amounts of ejective feedback that returns gas to 
the IGM before it forms stars. Such ejective feedback is also 
needed to explain the observed en richment of the IGM (e.g., 
lOppenheimer fc Pavel [200^ .1 20081 ). Ejected, metal-enriched 
gas will affect the cooling of the IGM in massive halos, and 
it will reduce the rate at which mergers add stars and gas 
to massive galaxies. We will investigate the impact of winds 
and metal-enrichment on galaxy growth in future work with 
simulations that include these effects. However, our present 
results suggest that the cooling of shock-heated, virialized 
gas, which has been the focus of many analytic models of 
galaxy growth spanning more than three decades, might be 
a relatively minor element of galaxy formation. 
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